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Abstract 

Background: Breeding livestock for improved resistance to disease is an increasingly important selection goal. 
However, the risk of pathogens adapting to livestock bred for improved disease resistance is difficult to quantify. 
Here, we explore the possibility of gastrointestinal worms adapting to sheep bred for low faecal worm egg count 
using computer simulation. Our model assumes sheep and worm genotypes interact at a single locus, such that 
the effect of an A allele in sheep is dependent on worm genotype, and the B allele in worms is favourable for 
parasitizing the A allele sheep but may increase mortality on pasture. We describe the requirements for adaptation 
and test if worm adaptation (1) is slowed by non-genetic features of worm infections and (2) can occur with little 
observable change in faecal worm egg count. 

Results: Adaptation in worms was found to be primarily influenced by overall worm fitness, viz. the balance 
between the advantage of the B allele during the parasitic stage in sheep and its disadvantage on pasture. Genetic 
variation at the interacting locus in worms could be from de novo or segregating mutations, but de novo mutations 
are rare and segregating mutations are likely constrained to have (near) neutral effects on worm fitness. Most other 
aspects of the worm infection we modelled did not affect the outcomes. However, the host-controlled mechanism 
to reduce faecal worm egg count by lowering worm fecundity reduced the selection pressure on worms to adapt 
compared to other mechanisms, such as increasing worm mortality. Temporal changes in worm egg count were 
unreliable for detecting adaptation, despite the steady environment assumed in the simulations. 

Conclusions: Adaptation of worms to sheep selected for low faecal worm egg count requires an allele segregating 
in worms that is favourable in animals with improved resistance but less favourable in other animals. Obtaining 
alleles with this specific property seems unlikely. With support from experimental data, we conclude that selection 
for low faecal worm egg count should be stable over a short time frame (e.g. 20 years). We are further exploring 
model outcomes with multiple loci and comparing outcomes to other control strategies. 



Background 

Sheep selected over many generations for low faecal 
worm egg count (WEC) show large reductions in WEC 
when infected with gastro-intestinal parasites [1-3]. This 
reduction in WEC indicates a reduced reproductive cap- 
acity, or reduced fitness, for worms in these hosts. A 
concern often voiced is that worm adaptation to sheep 
bred for low WEC may occur. The expectation is that 
natural selection results in worms increasing their fitness 
(or reproductive capacity) over time and this could erode 
the benefits of selecting sheep for low WEC. However, 
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to date these expectations have not been realised in 
experiments [4-7] and, anecdotally, breed differences in 
WEC appear to be stable across time. In this paper, we 
investigate whether the expectation of worm adaptation 
to low WEC sheep is realistic and attempt to determine 
factors that influence this potential rate of evolution. 

A compelling explanation of why worms are not 
observed to adapt to sheep selected for low WEC is lack- 
ing. Numerous authors have put forward possible explana- 
tions, including features of the host-parasite relationship 
that slow the rate of adaptation, such as the skewed distri- 
bution of WEC, which results in the majority of worms 
being harboured by susceptible hosts; and the insensitivity 
of WEC as a measure of adaptation e.g. [8-10]. However 
these assertions are rarely (if ever) quantified. It is unclear 
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which factors influence adaptation of worms to hosts that 
have been selected for improved resistance. More gene- 
rally, it is important to know which factors might also be 
applicable to other livestock disease systems (i.e. species 
and pathogens) such as selection of salmon or chickens 
for improved resistance to bacterial or viral infections. 
This investigation aims at identifying the key theoretical 
parameters that determine adaptation of nematodes in 
sheep. Results of this study may also give insight into 
the possible influence of these parameters in other 
situations. 

To explore factors that may affect the rate of parasite 
adaptation to host genotype, we describe a model that 
combines interacting descriptions of host genetics, para- 
site genetics and disease biology. We explore its key 
parameters and test hypotheses relating to observable 
rates of parasite adaptation. The model is parameterised 
using a single locus, gene-for-gene model. We inves- 
tigate the factors that generate or maintain genetic 
variation in worms at this interacting locus and test two 
specific hypotheses relating to worm adaptation i.e. (1) 
the relative rate of adaptation in worms is slowed by 
features of the host-parasite relationship, such as the 
skewed distribution of WEC and the different mecha- 
nisms by which sheep reduce WEC, and (2) WEC is an 
insensitive measure of worm adaptation. 

Methods 

Overview 

This paper describes results from two models, referred 
to as the 'full' and reduced' models. Both models have 
the same underlying interaction between genotype of the 
sheep and the worm but the full model also includes 
non-genetic factors related to the infection of sheep with 
gastrointestinal worms. Details of the genotype inter- 
action and the worm infection used in the full model are 
given in the section entitled "Model construction". Key 
definitions are that the A allele results in resistance' or 
low WEC in sheep, compared to the alternate a allele, 
and the B allele in worms is favourable for parasitizing 
sheep with the A allele (compared to the alternate b 
allele). In the section entitled "Model interrogation", we 
describe five scenarios, which are used to explore the 



values of the parameters that maintain genetic variation 
in worms (Scenarios 1 and 2), describe the properties of 
the worm infection and genotype by genotype inter- 
action (Scenarios 3 and 4) and test the hypotheses 
(Scenario 5). For simplicity, all simulations aim at 
maintaining a constant frequency of the A allele in sheep 
flocks over time and therefore we detect adaptation of 
worms to particular flocks by observing an increase in 
the frequency of the B allele in worms. To assess the im- 
pact of selecting sheep for lower WEC (i.e. an increase 
in the frequency of the A allele), we simulate a range of 
sheep flocks with different frequencies for A and deduce 
the impact on the worm population of selecting sheep 
for lower WEC by interpolation between the different 
flocks. 

Model construction 

Interaction between sheep and worm genotypes 

The model aims at describing worm adaptation to sheep 
with low WEC and hence it quantifies only differences 
between sheep genotypes for which worms possess the 
genetic potential to adapt. The interaction assumes bi- 
allelic loci with additive allelic effects (Table 1). The 
frequency of the A allele in sheep is x and the fre- 
quency of the B allele in worms is p. Genotype frequen- 
cies are given by Hardy- Weinberg proportions [11], 
such that sheep with AA genotype have a population 
frequency of x 2 . 

The B allele in worms is favourable for parasitizing 
sheep that carry the A allele (i.e. those animals with im- 
proved resistance to worms). The magnitude of this sur- 
vival advantage in worms is determined by the parameter 
s h (survival in the host). The model also has capacity for 
pleiotropic effects for the B allele so that the allele may be 
unfavourable for survival in free-living stages of the worm 
lifecycle (Table 1) [12]. The degree of pleiotropy is con- 
trolled by the free-living survival parameter s p (survival on 
pasture), such that the pleiotropic effects do not act when 
s p is zero. Therefore, there is the potential for p to increase 
when the worms are harboured by sheep with the A allele 
and for p to decrease during the free-living stage. There is 
a subtle but important distinction between worm fitness; 
defined as the success of an individual throughout its 



Table 1 Relative fitness of worm genotypes (per generation) as a function of host genotype 

Worm genotype Sheep genotype 



Susceptible {aa) Heterozygous {Aa) Resistant {AA) 





frequency 


(1-x) 2 


2x(1-x) 


x 2 


Wild-type {bb) 


(1-p) 2 


7 


1 - V4S h 


1 - V2S h 


Heterozygous {Bb) 


2p(l-p) 


1 - V2S P 


1 - V2S P 


1 - V2S P 


Alternate {BB) 


P 2 




(i+y 4 s h )(i-s p ) 


(i+y 2 s h )(i-s p ) 



s h = survival advantage of the alternate B allele for worms when harboured by sheep carrying the 'resistance' A allele; s p = survival trade-off (mortality) for the B 
allele, independent of sheep genotype. 
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lifecycle, and worm survival! which is relevant to only one 
aspect of the lifecycle. 

The reduced model 

The reduced model uses only population-level genetic 
information to predict the expected change in the fre- 
quency of the B allele in worms (per generation) for a 
given frequency of the A allele in sheep. The generational 
fitness (Q) is the average fitness of the worm bb, Bb and 
BB genotypes, weighted by the expected frequency of the 
three host genotypes from Table 1. That is: 



= (l-*) 2 + 2x(l-x) ( 1- ^s h 



l-\s h 



and develop over several days. Initially, newly hatched 
larvae from these eggs are immature but after 7 days the 
mature larvae are available for ingestion by the grazing 
flock. Thus, the model has subpopulations of worms 
defined by day age-classes, either on pasture or within a 
host, and a flock of grazing sheep. Each day age-class of 
worms has its own frequency of the B allele (p). 

We explore three postulated mechanisms for the sheep 
genotype to influence WEC in sheep: (i) reduce worm 
establishment, (ii) increase adult worm mortality and 
(iii) reduce adult egg production [3,13,14]. The A allele 
acts through these mechanisms in sheep, facilitating the 
potential for worms to adapt in the full model. For each 
mechanism, the A allele reduces WEC by a defined 
amount a (described below), such that homozygous AA 
sheep have a WEC that is 2a lower than the WEC of 
homozygous aa sheep. In the simulations, each mecha- 
nism to reduce WEC is tested independently. 



(l-s v 



(1-x) 2 + 2x(l-x) ( 1 + -s h 
4 



The change in worm allele frequency [11], is then: 

P 2 Q B b + p(l-p)Q.Bb \ 



P 2 ei BB + 2p(i-p)ciBb + (i-^) 2 U 



(1) 



bb j 



where p and p' are the frequency of the B allele in the 
current and next worm generations. In contrast to the 
full model, the reduced model assumes an infinite popu- 
lation of worms and uses discrete generations to make 
deterministic predictions of the change in p. 

The full model 

The full model includes a flock of sheep grazing on a 
common pasture. Each sheep has its own level of resist- 
ance and worm population, with worms continuously 
infecting and re-infecting sheep. The worm lifecycle 
obeys the time-dependent dynamics seen with worm 
infections and, over a season, animals slowly develop 
increased immunity to the worm infections. The flock 
has a skewed distribution of WEC measurements, in 
accordance with field observations. 

The full model is described in full detail in Additional 
file 1. Briefly, the sheep by worm genotype interaction is 
expanded into a series of daily steps. Each day, sheep 
ingest worm larvae and initiate a set of interactions 
unique to the infected animals genotype and the geno- 
type of the larvae ingested on that day. Over time, a 
proportion of the ingested larvae develop into adult 
worms. Adult (female) worms produce eggs until they 
eventually die. The eggs produced by worms on a given 
day can be observed in each animal via a measurement 
of WEC. Eggs are deposited on the pasture in the faeces 



Model interrogation 

Genotype by genotype interaction and the maintenance of 
genetic variation at the interacting locus in worms 

Worm alleles at the interacting locus could either be de 
novo mutations or variants that already segregate in the 
population. The reduced model was used to examine 
these two situations and the properties of the genotype 
by genotype interaction. There are two scenarios, as 
defined in Table 2, which use a range of parameter 
values for p 0 (the initial frequency of the B allele in 
worms), s h (the survival advantage of worms with the B 
allele in sheep with the A allele), s p (the survival penalty 
on pasture for the B allele) and x (the frequency of the 
A allele in sheep). Scenario 1 investigates the impact of 
changes in p 0 on the outcome for p over many genera- 
tions, while other parameters are held constant. This 
scenario is relevant to study the properties of de novo 
mutations that are initially very rare in the population. 
Scenario 2 explores the maintenance of segregating vari- 
ants and the impact of changing the survival advantage 
(s h ) and trade-off (s p ) of the B allele in worms. The 
frequency of the A allele in sheep was equal to 0.2, 0.5 
or 0.8. Interpolation across simulations in Scenario 2 
shows the properties of the segregating B allele for 
worms in flocks of sheep with an increasing frequency 
of the A allele. 

Evaluation of the worm infection under the full model 

The properties of the worm infection for the full model 
were examined in two scenarios under conditions 
whereby the sheep by worm genotype interaction could 
not influence the results. These scenarios aimed at 
demonstrating the properties and appropriateness of the 
simulated worm infection. Parameters for these scena- 
rios are defined in Table 2 (Scenarios 3 and 4). The aim 



Kemper et al. Genetics Selection Evolution 2013, 45:14 
http://www.gsejournal.Org/content/45/1 /1 4 



Page 4 of 14 



Table 2 Summary of the scenarios tested 



Scenario Number of worm generations s h (%) 
or years 8 


s p (%) 


Po 


X 


Model type b 


Mechanism to reduce 
sheep WEC C 


Evaluation of the genotype-by-genotype interaction 












1 80 gen 20 


0 


0.001, 0.01, 
0.1,0.3 


1.0 


reduced 


- 


2 80 gen 0 to 30 


Oto 30 


0.3 


0.2, 0.5, 0.8 


reduced 




Evaluation of the host-parasite interaction properties 












3 1 yr 




0 


0 


full 




4 1 yr 25 


10 


0.5 d 


0.5 


full 


E, M, F 


Testing hypotheses 7 and 2 












5 20 yr 25 


10 


0.3 


0.0 to 1 .0 


reduced, full 


E, M, F 



a the reduced model is defined by worm generations and the full model is defined by years of simulation; results show that there are 2.2 worm generations per 
year of simulation; s h = survival advantage of the B allele in sheep; s p = survival trade-off on pasture for the B allele; p 0 = initial frequency of the B allele in worms; 
x = frequency of the low WEC A allele in sheep; b see text for details of the full model, which includes complexities of the host-parasite relationship, and the 
reduced model, which explicitly excludes these factors; c postulated mechanisms to reduce WEC in sheep through the 'resistance' A allele: E = reduced worm 
establishment, M = increased worm mortality, and F = reduced female worm fecundity; d to evaluate the effect of the low WEC allele for each mechanism of the 
sheep 'resistance' A allele, the worm allele frequency was fixed (p = 0.5) during the transition stages of the worm lifecycle. 



of Scenario 3 was to evaluate the characteristics of the 
worm infection in a single season. This scenario used a 
susceptible' sheep flock (i.e. x = 0) and p 0 was zero, 
therefore there was no potential for the worm popula- 
tion to adapt. Scenario 4 tested the properties of the A 
allele in sheep to ensure that each mechanism to reduce 
WEC had an equivalent effect. That is, we wanted to 
calibrate the mechanisms such that homozygous AA 
sheep had a WEC that was 2a lower than homozygous 
aa (susceptible') sheep. Special conditions were required 
for this scenario to prevent accumulation of effects 
throughout a season, in order to ensure equivalent test 
conditions between the mechanisms at the end of the 
season, and also because the effect a of the A allele 
depends on the frequency of the B allele in worms. 
Effects can accumulate within a season because the 
genotype profile of the worm population changes continu- 
ally, affecting subsequent WEC measurements. Thus, the 
allele frequency of the worm population was fixed to 0.5 
during each transition stage (i.e. when ingested as L 3 
larvae in sheep and when first deposited as eggs on 
pasture). A frequency of 0.5 was chosen to ensure equal 
numbers of homozygotes for worm genotypes. The ave- 
rage allele effect a was calculated as 

where WECaa is the average WEC of sheep that are 
homozygous for the low WEC allele, and WEC aa is that of 
homozygous susceptible sheep. 

Testing hypothesis 1: non-genetic factors in the worm 
infection cycle can slow the rate of worm adaptation 

Hypothesis 1 proposes that features of the host-parasite 
interaction will slow the rate of adaptation. These factors 



are either (a) factors involved in worm infection (e.g. a 
skewed distribution of WEC or overlapping worm gene- 
rations) or (b) the mechanisms of resistance. This hy- 
pothesis was tested by assessing if the reduced model 
could predict the change in p from the full model, using 
each of the three postulated mechanisms to reduce 
WEC. Thus, we tested all factors included in the full 
model that are not involved in the genotype by genotype 
interaction. Parameters to test the hypothesis are 
detailed in Scenario 5 (Table 2). The parameters were 
defined such that worms with the B allele can have a 
25% survival in sheep that carry the A allele but there is 
a 10% survival trade-off for the B allele on pasture (i.e. 
s h = 0.25, s p = 0.1). Sheep with one copy of the A allele 
had approximately a 25% reduction in WEC. Each 
proposed mechanism by which sheep reduce WEC was 
tested in turn, for scenarios in which the frequency of 
the A allele was increased incrementally from 0 (suscep- 
tible) to 1 (resistant). The initial frequency of the B allele 
in worms (p 0 ) was 0.3. Results presented for these 
scenarios were the mean of five replicate sheep flocks. 

Testing hypothesis 2: WEC is an insensitive measure of 
worm adaptation 

Worm adaptation to sheep with the low WEC A allele 
may be expected to result in an increase in WEC over 
time for these animals, i.e. when worms adapt to sheep 
with the A allele, they increase their egg output and 
hence WEC. In our model, we expect an increase in p 
(i.e. adaptation') to be associated with an increase in 
WEC. We used WEC as an indicator of adaptation for 
data generated from Scenario 5. A linear model (y = a + bt) 
was fitted to the replicate mean WEC over time. Adap- 
tation to the A allele in sheep was declared if the 
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regression coefficient (b) was significantly greater than 
zero (P < 0.001). 

Results 

Model parameters for the sheep-by-worm genotype 
interaction and the maintenance of genetic variation at 
the interacting locus in worms 

The reduced model was used to investigate the effect of 
changes to the underlying parameters for the genotype 
interaction; viz. changes to p 0 , s p , s h and x (i.e. the initial 
frequency of the B allele in worms, the survival advan- 
tage and trade-off for the B allele and the frequency of 
the "resistance" A allele in sheep). The results from 
Scenario 1 are important when considering de novo 
mutations as a source of genetic variation for adaptation 
in worm populations. This is because de novo mutations 
are expected to be very rare alleles (i.e. a very low p 0 ). 
The results indicate that changing p 0 had little impact 
on the final outcome for the B allele but that p 0 did 
affect the number of generations required before an 
observable change in p occurred (Figure 1). Thus, 
worms adapt (p tends to 1) when the B allele increases 
worm survival in sheep with the A allele and when the B 
allele has no pleiotropic effects, regardless of the initial 
frequency of the B allele. Changing p 0 affected the total 
number of generations required until fixation but had 
no affect on the rate of change at a given allele fre- 
quency. For example, an allele that was initially very rare 
(p 0 = 0.001) took 80 generations to achieve fixation 
compared to an allele that was initially at a higher fre- 
quency (i.e. p 0 = 0.1). Thus, de novo mutations are 
expected to behave like other alleles with similar pro- 
perties (i.e. Sh and s p ) but observable changes for these 
alleles will be slower compared to similar alleles at 
higher frequencies. 

c \ 



1.0 - 



Q. 








nc 


0.8 - 


CD 




=3 








fre 


0.6 - 


CD 




0) 




"ro 


0.4 - 


3 








c 


0.2 - 


£ 










0.0 - 



0 20 40 60 80 

generations of worms 

Figure 1 Changes in the frequency of the alternate B allele in 

worms when the initial frequency of B is either (i) 0.001 7 (ii) 

0.01 , (iii) 0.1 or (iv) 0. Under Scenario 1, there is a 20% survival 

advantage {s h ) for the B allele in sheep carrying the A allele and no 

mortality trade-off on pasture (s p = 0). 
\ ) 



Scenario 2 investigated changing the magnitude of the 
survival effects (s h and s p ) of the B allele. Changes in 
these parameters had dramatic impacts on the final out- 
come for the B allele. Figure 2 summarises results from 
many simulations using the reduced model, when the A 
allele in sheep had a frequency of 0.2, 0.5 or 0.8. The 
results indicate that for most of the tested parameter 
combinations, worms adapt either by fixing (p = 1) or 
completely losing (p = 0) the B allele. If the survival 
trade-off of the B allele on pasture is greater than its 
survival advantage in sheep, then the B allele will be lost. 
Since x < 1, worm survival on pasture is always more 
important than worm survival in the host and thus in 
most cases when s p > s h , then p tends to 0. This is 
because the B allele is advantageous to worms only when 
they are harboured by sheep with the low WEC A allele 
(which is uncommon if the frequency of A is low), while 
the survival trade-off for the B allele acts on all larvae 
while on pasture. Thus, the relative importance of 
survival on pasture depends on the frequency of the A 
allele in sheep. For example, a mutation in worms with a 
30% survival advantage in sheep (s h = 0.3) and a 15% 
survival trade-off on pasture (s p = 0.15) would be lost when 
the frequency of the A allele in sheep was low (x = 0.2), but 
would remain segregating if the frequency of the A allele 
was higher (i.e. x = 0.5). 

If the mutations available for adaptation come from 
segregating variants at the interacting locus in worms, 
they must persist after many generations of natural 
selection. Results from Figure 2 show the type of vari- 
ants that segregated in the worm population after 80 
generations of natural selection. These results indicate a 
grey area of segregating variants where the B allele in 
worms is neither fixed nor lost. However, given (e.g.) 
200 or 1000 generations, this grey area of segregating 
alleles would be even smaller. For polymorphism to be 
maintained indefinitely, the homozygous BB and bb 
worm genotypes are required to have equal fitness (i.e. 
p' = p from equation (1) only when s p = l- Qbb / n ) . 
Thus over many generations, only alleles where 
the survival advantage is precisely balanced by the 
survival trade-off remain segregating. Some unbalanced 
alleles may also appear to segregate (or drift) because as 
Sp + abb / DiBB tends to 1, the rate of change in the fre- 
quency of the B allele tends to 0. If adaptation in worms 
occurs by a change in frequency of segregating variants, 
then these alleles must affect worm survival but must 
also have a neutral or near neutral overall effect on 
worm fitness. Otherwise natural selection would act to 
fix the favourable allele in worms regardless of the worm 
resistance exhibited by the sheep. The values for the 
survival advantage and the trade-off for alleles which 
segregate indefinitely depend on the frequency of the 
low WEC A allele in sheep. As the frequency of the 
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€ A B C 

0 0.1 0.2 0.3 0 0.1 0.2 0.3 0 0.1 0.2 0.3 

mortality on pasture, s p 

Figure 2 Frequency of the alternate B allele in worms after 80 generations when the interacting A allele in sheep is (A) 0.2, (B) 0.5 or (C) 

0.8. Under Scenario 2, the B allele survival advantage in sheep carrying the A allele (s hl y-axis) and pasture mortality (s p , x-axis) range from 0 to 30%. 



A allele increases, the relative importance of survival on 
pasture decreases (i.e. the polymorphic combinations in 
Figure 2 rotate clockwise from vertical, s p ~ 0, to 45° 
from vertical, s p tends to s h , as the frequency of A in- 
creases). Interestingly, the number of parameter sets for 
which alleles remain segregating (i.e. the size of the grey 
area) does not depend on the frequency of the A allele. 
For example, the grey area is an equal size when only 
one (x = 1.0) or all three (e.g. x = 0.2) sheep genotypes 
are present, implying that heterogeneity in sheep does 
not increase worm heterogeneity in our model. 

The B allele at the interacting locus only remains seg- 
regating after many generations when there is no overall 
fitness effect for the allele; either the survival advantage 
of the B allele in sheep with the A allele is balanced by 
the survival trade-off on pasture, or the B allele does not 
affect worm survival. However, this balance depends on 
the frequency of the A allele. Thus, if we take a given 
sheep population and imagine selecting for an increase 
in frequency of the A allele (i.e. increased resistance to 
worms or a reduced WEC), then worm alleles that previ- 
ously had no fitness effect may become favourable and 
begin to increase in frequency. The rate of increase in 
the B allele is determined by the magnitude of change in 
the frequency of the A allele, i.e. the change in the 
frequency of the A allele determines how unbalanced 
the survival advantage and trade-off in worms becomes. 
If the frequency of the A allele is slightly increased, then 
there is only a small change in the selection pressures 
that act on the worms (i.e. a small rotation in the grey 
area) and, over time, the frequency of the B allele 
increases slowly. If there is a large increase in the 
frequency of the A allele, then the change in the selec- 
tion pressures for the worm population is greater and 
this will cause a more rapid increase in the frequency of 



the B allele. In summary, if adaptation in worms is to be 
observed from segregating variants, the mutation needs 
to change from having a neutral to having a favourable 
effect on worm fitness. The magnitude of this change in 
fitness determines the rate of adaptation, i.e. the rate of 
increase in frequency of the favourable B allele. 

Properties of the worm infection from the full model 

Scenario 3 aims at investigating the epidemiological 
properties of the full model without the complication of 
genetic variation in either the sheep or worm popula- 
tions. The degree of pasture contamination shows the 
presence of a pasture reservoir of worm larvae (Figure 3). 
There are two annual peaks in total pasture contami- 
nation, the first at around day 45 and the second at day 
85; and two slightly earlier peaks in new infective L 3 
stage larvae on pasture. These peaks are caused by two 
successive generations of worms infecting the sheep. 
WEC and worm burden (not shown) show skewed dis- 
tributions, which is characteristic of these traits. Mean 
adult burden was 19 310 and mean WEC 769 eggs per 
gram (epg). The degree of skewness in WEC agrees with 
previous experimental results [7]. The average gene- 
ration length for the worm population was 45 days. 
Hence one year (or 100 days of simulation) is equivalent 
to 2.2 worm generations, in accordance with field obser- 
vations [15,16]. These results indicate that our model 
adequately captured the key epidemiological characteris- 
tics of worm infections that we aimed to include in our 
model. 

The effect of the low WEC A allele in sheep and the 
corresponding change in the frequency of the B allele in 
worms is reported in Table 3, taking a census at day 100. 
For the purpose of testing the model, we held the fre- 
quency of the B allele for ingested larvae and new eggs on 
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pasture at 0.5, preventing accumulation of changes to 
worm allele frequency within a season. This procedure was 
only used when testing the model to ensure correct para- 
meterisation. Depending on the mechanism implemented, 
sheep with one copy of the A allele either reduce the pro- 
portion of larvae that establish by 0.029 (untransformed 
scale), increase daily mortality by 0.038, or reduce female 
fecundity by 16 eggs per worm. The low WEC allele (A) 
reduced WEC at day 100 by approximately 25% (or 100- 
150 epg) for all mechanisms; hence the WEC of homozy- 
gous AA sheep was 40 to 60% of the WEC of homozygous 
aa sheep. Table 3 demonstrates some important conse- 
quences of the mechanisms to reduce WEC and of the 
epidemiological model. These include: 



A reduction in the number of adult worms as x 
increases (from aa to Aa to AA sheep), when the A 
allele reduces worm establishment or increases adult 
mortality. The adult worms are younger when 
harboured by sheep with an A allele that acts by 
increasing adult worm mortality. As expected, there 
is no change in the number of adult worms in sheep 
with the A allele when it acts to reduce worm 
fecundity. 

There is an increase in the frequency of the B allele 
for adult worms harboured by sheep with the A 
allele when the mechanism of the A allele is to 
reduce worm establishment or increase adult worm 
mortality. Similar to point 1 above, no differences in 



Table 3 Effect of the low WEC A allele in sheep depending on the mechanism by which sheep reduce WEC 



Mechanism of 
the A allele 



Sheep 
genotype 



Mean adult worms 3 



WEC mean b 



Count 



Age (days) Count (epg) 



Effect of 
allele A (a%) 



Mean pasture L 3 



Max nb (xlO 7 ) 



Age (days) 



E 


aa 


11 769 


0.500 


31.3 


487 


0.500 










Aa 


8 776 


0.516 


31.2 


352 


0.516 


22.4 


4.72 


0.489 15.2 




AA 


6610 


0.531 


31.2 


268 


0.531 








M 


aa 


12 581 


0.500 


30.7 


501 


0.500 










Aa 


8 916 


0.518 


27.3 


356 


0.518 


29.7 


4.93 


0.488 15.2 




AA 


5 067 


0.534 


25.6 


203 


0.534 








F 


aa 


11 134 


0.500 


31.3 


449 


0.500 










Aa 


11 344 


0.500 


31.3 


354 


0.508 


21.9 


4.22 


0.488 15.2 




AA 


11 267 


0.500 


31.3 


252 


0.516 









Postulated mechanisms to reduce WEC in sheep through the 'resistance' A allele: E = reduced worm establishment, M = increased worm mortality and F = 
reduced female worm fecundity; epg = eggs per gram; a number, age, and frequency of the alternate B allele (p) for adult worms; b average effect on WEC (a) and 
p for deposited eggs; c number, frequency of p and age for infective L 3 larvae on pasture. 
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the frequency of the B allele are evident between 
AA, Aa and aa sheep genotypes when the A allele 
acts by reducing worm fecundity. 

3. The frequency of the B allele is unchanged from 
adult worms to eggs when the A allele acts by either 
reducing worm establishment or increasing adult 
mortality. However, the frequency of the B allele is 
greater in the eggs than in adult worms when the A 
allele acts by reducing worm fecundity. 

4. The frequency of the B allele increases in sheep 
(i.e. in sheep with the A allele) but decreases during 
the pasture stage of the lifecycle. That is, there is a 
survival trade-off on pasture for the B allele. 

5. The magnitude of increase in the frequency of the B 
allele differs between the three mechanisms of the A 
allele. When the A allele acts by reducing adult 
worm fecundity, the change in allele frequency is 
approximately half of that when the A allele acts by 
reducing establishment or increasing adult worm 
mortality. This is insightful for the fecundity 
mechanism. When the A allele acts by reducing 
worm fecundity, it affects only the female worms 
that contribute to the next generation and not the 
male worms (i.e. the effect is sex-limited). Thus, the 
magnitude of increase for p in the eggs is only half 
of that when the mechanism of the A allele acts on 
both sexes. 

Testing hypothesis 7: non-genetic features of the host- 
parasite relationship slow the rate of worm adaptation 

The change in frequency of the worm B allele is shown 
over 20 generations for Scenario 5 in Figure 4 (plots A, 
B, C). For most frequencies of the A allele in sheep, the 
B allele in worms either moves towards fixation or loss. 
The rate of change (i.e. the slope of the line) depends on 
the frequency of the A allele in sheep. Thus, as the sheep 
population has lower WEC (i.e. as x increases), the slope 
describing the change in the frequency of the B allele in 
worms is increasingly positive, indicating that the B 
allele is favourable. Adaptation was most rapid when the 
low WEC A allele in sheep was fixed (x = 1.0). 

Comparing the change in allele frequency in worms 
between the outputs of the four models allows the first 
hypothesis to be tested (Figure 4). If non-genetic factors 
involved with the worm lifecycle, such as the skewness 
of WEC and a reservoir of infective larvae on pasture, 
had an effect on the rate of worm adaptation, then there 
would be observable differences between outputs from 
the full and the reduced models. Second, if differences in 
the rate of adaptation are due to the mechanisms by 
which WEC is reduced in sheep, there would be diffe- 
rences between the three sets of results from the full 
model. Figure 4 shows a strong similarity between the 
allele frequency change for the establishment and 



mortality mechanisms from the full model and the re- 
duced model, bearing in mind that one year approxi- 
mates 2.2 parasite generations. There may be a slightly 
higher final frequency of the B allele for worms that 
parasitize sheep with the mortality mechanism, com- 
pared to results from the establishment mechanism or 
from the reduced model. This is probably due to the 
non-linearity of the mortality mechanism and the accu- 
mulation of effects over time due to changes in the age 
structure of the worm population. In contrast, the re- 
sults from the reduced model do not predict the results 
from the full model when the mechanism to reduce 
WEC is by decreasing adult worm fecundity. For this 
mechanism, the rate of change in the frequency of the B 
allele in worms was negative at most frequencies of the 
A allele in sheep. Hence the B allele is mostly unfavour- 
able for the worms. This arises because selection pres- 
sure is reduced for this mechanism, as observed in 
Table 3. The results from Figure 4 suggest that most fea- 
tures of the host-parasite interaction that were consid- 
ered in the full model, including the between-host 
distribution of worms, have little impact on the rate of 
adaptation in worms. However, when sheep reduce 
WEC by reducing worm fecundity, a slower rate of 
adaptation in worms is predicted. 

Notably, although the rate of change in the frequency 
of the B allele in worms is lower when WEC is reduced 
through the fecundity mechanism compared to other 
mechanisms, the rate of loss for the B allele in the 
susceptible population (x = 0) is similar for all four sets 
of results. This property is best explained with reference 
to the gene-by-gene interaction defined in Table 1. If 
only susceptible aa sheep are considered, then worm 
fitness is only determined by the value of s p . This para- 
meter is constant for this scenario regardless of the 
mechanism and hence Figure 4 shows the same rate of 
loss for all unselected sheep. The sex-limited action of 
the fecundity mechanism only affects the relative fitness 
of the worm genotypes when x is greater than zero. In 
effect, the sex-limited action halves the values of s h in 
Table 1. This changes the balance for fitness between 
worm genotypes, leading to results for the fecundity 
mechanism being different from results for the other 
mechanisms. 

Testing hypothesis 2: WEC is an insensitive measure of 
worm adaptation 

The flock average WEC measurements showed the clear 
effect of a reduction in WEC with increasing frequency 
of the A allele in the sheep population for each of the 
resistance mechanisms (Figure 5 and Table 4). However, 
the WEC for sheep flocks in which the A allele was fixed 
(i.e. x = 1.0) was much lower than that estimated by the 
effect of the resistant allele. That is, the homozygous AA 
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Figure 4 Mean frequency of the B allele in worms during 20 years of simulation using either the full (A-C) or reduced (D) model. Under 
scenario 5, there is a 25% survival advantage {s h ) and 10% pasture mortality (s p ) for the B worm allele. The mechanisms by which sheep reduce 
faecal worm egg count (WEC) are (A) reduced worm establishment, (B) increased worm mortality, or (C) reduced worm fecundity; the frequency 
of the low WEC A allele in sheep is investigated at 1 1 levels, from low (light blue) to high (dark blue) frequencies and also when the allele is 
absent (black). 



genotype sheep had a WEC of 40% to 60% of the aa 
sheep in Table 3 but the WEC of the AA flock was 20 to 
30% of the aa flock (using the average of years 1 to 5) in 
Scenario 5. This difference is primarily driven by the 
within-year effects on pasture contamination. For ex- 
ample, the susceptible' sheep flocks from Table 4 have 
an initial WEC of -700 epg, while the aa animals from 
Table 3 have a WEC of -475 epg. The susceptible flock 
from Scenario 5 has a higher WEC because pasture con- 
tamination for the flock is higher (i.e. there are no sheep 



with the AA genotype, while in Scenario 4 all three 
genotypes contribute to pasture contamination) and the 
mortality on pasture for worms is lower because the fre- 
quency of the B allele in worms is lower (i.e. p 0 is 0.3 in 
Scenario 5 compared to 0.5 in Scenario 4). In Scenario 
5, an additional - 25% reduction in WEC is observed in 
an AA flock (x = 1) because of reduced worm burden 
and larval intake arising from lower larval 
contamination from low WEC sheep and increased mor- 
tality for B allele worms on pasture. 
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Figure 5 Mean faecal worm egg count (WEC) during 20 years of simulation of worm infection in sheep under Scenario 5. Mechanism by 
which sheep reduce WEC are (A) reduced worm establishment or (B) increased worm mortality; the frequency of the low WEC A allele in sheep 
was investigated at 1 1 levels, from low (light blue) to high (dark blue) frequencies and also when the allele is absent (black). 



Notably, the WEC in flocks for which obvious adapta- 
tion occurred (e.g. the x = 1.0 flock from Figures 4B and 
5B) did not recover to the level seen in aa sheep flocks 
(x = 0). This is because the B worm allele suffers a 
survival penalty on pasture in our simulations, lowering 
the number of infective larvae on pasture when the 
worm population has an increased frequency of the B 
allele. Thus, although worms improve survival by adap- 
ting to sheep with AA genotype, overall fitness for 
worms with the B allele remains lower than wild-type 
worms in susceptible sheep because of the allele s pleio- 
tropic effects. 

Changes in WEC over time were generally a poor indi- 
cator of adaptation because of its high annual fluctua- 
tions (Figure 5). The annual fluctuations are not due to 
variation in environmental factors such as climate (these 
parameters were constant) but are likely caused by 
random processes in the model, such as daily variation 
in feed intake. Although the significant regression of 
WEC on time indicated quantifiable adaptation, reflec- 
ting the increase in the frequency of the B allele 
(Table 4), the changes in WEC were rather small and 
noisy. For example, significant adaptation was indicated 
when the A allele reduced establishment for the x = 0.1 
and 0.3 flocks but not for the intermediate x = 0.2 flock. 
As above, the increase in WEC was sometimes asso- 
ciated with a loss of the B allele because of its pleio- 
tropic effects on pasture survival. It is unlikely that this 
magnitude of change in WEC would be observable 



under field conditions, where the level of larval chal- 
lenge varies across years, and changes are further 
masked by climatic fluctuations. 

Discussion 

This paper examined the conditions required for adapta- 
tion of worms at a single locus to sheep with an allele 
that confers low WEC, and tested two hypotheses 
related to rates of worm adaptation. For the former, our 
model was set up with specific conditions to facilitate 
adaptation of worms to host genotype by including 
interactions between sheep and worm genotypes (i.e. 
Table 1). This is necessary because when an allele in 
the worm consistently has favourable (or unfavourable) 
effects and no interaction with the genotype of the 
sheep, it would always be under positive (or negative) se- 
lection and, thus, such an allele could not be responsible 
for adaptation to specific sheep genotypes (i.e. low WEC 
sheep). Hence, we used a gene-for-gene model where 
the genotypes of sheep and worms could interact. The 
second requirement for adaptation is that worms need 
genetic variation at this locus. This variation could be 
through either de novo or segregating variants. However, 
de novo mutations are, by definition, initially rare alleles 
and we found that they may require many generations in 
the worm population before an appreciable change in 
allele frequency is observable. Furthermore, these alleles 
are quickly subjected to natural selection and purged 
if unfavourable. We found that the mutations that 
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Table 4 The regression of replicate mean WEC on time 
for flocks in Scenario 5 



Mechanism Frequency Flock average Regression 



Ap 



f the 
allele 


of the A 
allele (x) 


WEC (epg) 


coefficient 
(epg/yr) 








Year 
1 


Year 
20 


Estimate 


s.e. 




E 


0.0 


710 


710 


0.72 


0.87 




E 


0.1 


640 


670 


3.21 


0.80 


Decrease 


E 


0.2 


540 


510 


-1.07 


0.58 




E 


0.3 


470 


500 


2.08 


0.50 


Decrease 


E 


0.4 


380 


430 


1.16 


0.45 




E 


0.5 


340 


340 


-0.17 


0.50 




E 


0.6 


290 


300 


0.39 


0.36 




E 


0.7 


250 


260 


1.02 


0.37 




E 


0.8 


210 


230 


1.62 


0.26 


Increase 


E 


0.9 


180 


220 


1.83 


0.18 


Increase 


E 


1.0 


140 


190 


2.54 


0.20 


Increase 


M 


0.0 


720 


800 


2.94 


0.56 


Decrease 


M 


0.1 


640 


650 


1.22 


0.65 




M 


0.2 


580 


620 


2.39 


0.66 




M 


0.3 


500 


480 


-0.08 


0.62 




M 


0.4 


440 


440 


-0.89 


0.54 




M 


0.5 


370 


370 


-1.40 


0.43 




M 


0.6 


340 


350 


1.81 


0.39 


Increase 


M 


0.7 


290 


310 


1.07 


0.27 


Increase 


M 


0.8 


260 


320 


3.62 


0.24 


Increase 


M 


0.9 


220 


300 


4.35 


0.30 


Increase 


M 


1.0 


180 


290 


5.53 


0.29 


Increase 


F 


0.0 


720 


800 


2.94 


0.56 


Decrease 


F 


0.1 


660 


660 


0.10 


0.80 






0.2 


570 


540 


-1.13 


0.58 






0.3 


480 


520 


4.02 


0.63 


Decrease 




0.4 


410 


430 


-0.52 


0.53 






0.5 


360 


340 


0.08 


0.53 






0.6 


310 


320 


-0.69 


0.47 






0.7 


250 


240 


-0.33 


0.45 






0.8 


230 


210 


0.55 


0.37 






0.9 


180 


170 


-0.30 


0.29 






1.0 


140 


170 


0.99 


0.16 


Increase 



Postulated mechanisms to reduce WEC in sheep through the 'resistance' A allele: 
E = reduced worm establishment, M = increased worm mortality and 
F = reduced female worm fecundity; epg = eggs per gram; s h = 25%; s p = 10%; 
bold characters = coefficients significantly different to zero (P < 0.001) with the 
directional change in frequency for the alternate B allele in worms (Ap). 



remained segregating after several generations of expo- 
sure to sheep had very specific values for their survival 
advantage and trade-off. Only mutations for which the 
survival advantage in sheep was precisely balanced by 



their detrimental pleiotropic effects remained segrega- 
ting in the long term. Our model assumed the detrimen- 
tal effects comprised reduced survival on pasture but, 
equally, these effects could have occurred elsewhere in 
the lifecycle of the worm. Thus the basic requirements 
for adaptation in our model were presence of loci in the 
host and worm that interact and genetic variation in 
both the worm and host locus. 

The presence of genetic variation suitable for adap- 
tation of worms to sheep with improved resistance is 
unknown. A high degree of molecular polymorphism is 
observed in worm populations [17,18] and this variation 
could be equated with adaptive potential. Our results, 
particularly from Scenario 2, suggest that polymor- 
phisms in worms with appreciable frequencies (and 
therefore old' mutations segregating in the population) 
are likely neutral, or near neutral, for overall worm fit- 
ness. This is because worms have infected sheep for 
many generations and natural selection will have had 
opportunity to maximise worm fitness. Polymorphisms 
that segregate in worm populations could affect fitness if 
they were younger, and thus exposed to fewer genera- 
tions of natural selection, or if, for example, they have 
transient effects on fitness [19]. However, the persistence 
of these polymorphisms will be inversely related to their 
net effect on fitness because any allele with strong bene- 
ficial effects will be rapidly fixed. Our model explored 
antagonistic pleiotropy as a mechanism to maintain 
polymorphism in worms and, although experimental evi- 
dence for antagonistic pleiotropy at individual loci is 
scarce [20], Jorgensen [12] reports reduced survival 
(larval development) for worms at the population level 
in a line of sheep selected for low WEC. The presence of 
polymorphisms in worms with antagonistic effects on 
survival is unknown but our model demonstrates that 
such a mechanism may be required to maintain poly- 
morphism at a locus with properties suitable for adapta- 
tion of worms to sheep genotype. Without pleiotropy, 
natural selection rapidly purges unfavourable alleles 
from the worm population and genetic variation is lost 
at the locus. We expect that most of the observed 
molecular polymorphisms in worms will have small 
effects on worm fitness but recognise that segregation of 
alleles with antagonistic pleiotropy or of young muta- 
tions with near-neutral effects on fitness remain a 
possibility. 

The final requirement for worm adaptation to sheep 
genotype is selection pressure. That is, there must be a 
difference in survival between worm genotypes that 
depends on sheep genotype. In our model, selection 
pressure depended on the survival advantage and trade- 
off for the allele at the interacting locus in the worm 
(allele B) and the allele frequency at the interacting locus 
in sheep (allele A). Our model behaved in predictable 
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and consistent ways across all scenarios with respect to 
worm adaptation (i.e. the change in the frequency of the 
interacting worm allele, p). Results from Scenario 5 
showed that, for a specific mutation, flocks with extreme 
frequencies of the interacting allele in the host (x = 0 or 
x = 1 for allele A) led to the most rapid rate of adapta- 
tion. This was because, under the conditions of Scenario 
5, the B allele remained segregating in the worm popula- 
tion after several generations when x ~ 0.5 and frequen- 
cies of either x = 0 or x = 1 caused the most disruption 
in this balance. The requirements for an interacting 
locus, i.e. genetic variation and selection pressure, 
should also apply to other host-pathogen systems. 

Given the 'ideal' starting conditions for worms to 
adapt, we were able to test two specific hypotheses with 
our model. The first hypothesis tested if the (non-gen- 
etic) features of the host-parasite interaction slowed the 
observed rate of adaptation. Some features, including a 
skewed distribution of WEC in sheep, overlapping gen- 
erations of worms and a reservoir of worm larvae on 
pasture, primarily impacted the generation interval of 
worms. Thus, after calibrating the outputs by generation 
interval in the worm, a simpler reduced model was able 
to accurately predict the results from the full model that 
included key host-parasite interactions. The reduced 
model showed that the outcome for the B allele in 
worms is essentially determined by the selection pres- 
sure per generation, i.e. the selection differential for each 
sheep genotype and the proportion of worms exposed to 
each sheep genotype, and by the generation length for 
the worm. Accurate estimates of these two parameters 
are important to apply these findings to other host- 
pathogen systems, although in practice this is not a triv- 
ial exercise. 

An important outcome of the model was that limiting 
female fecundity as the mechanism by which selection 
of sheep reduces WEC was found to halve the selection 
pressure on the worm population compared to other 
mechanisms for WEC reduction. This was because of 
the sex-limited action of this mechanism, whereby a 
female worm that is able to overcome the A allele in 
sheep has its favourable genotype diluted by mating with 
unselected male worms. Experimental evidence also 
suggests an important role for reduced worm fecundity 
in low WEC sheep [3,21-24], where WEC in sheep with 
immature immune responses seems to be suppressed 
by reduced worm fecundity rather than a (partial) reduc- 
tion in total number of worms. Other host-parasite 
systems also suggest slower adaptation when hosts 
reduce pathogen reproduction (or infection rate) as 
compared to pathogen numbers. For example, modelling 
in malaria shows imperfect vaccines that limit pathogen 
growth (i.e. numbers) are more likely to cause an 
evolutionary response in the pathogen compared to an 



infection-blocking mode of action [25]. Quantitative 
support for this prediction is evident in the successive 
adaptation of the Marek's disease virus to repeated re- 
leases of vaccines that limit tumor growth rather than 
viral reproduction or propagation [26,27]. Here, we 
make a similar favourable conclusion about mechanisms 
that limit pathogen reproduction, as compared to pa- 
thogen numbers, but our results are obtained through 
the sex-limited action of the fecundity mechanism, in 
contrast to other studies, which exploited within-host 
competition between pathogen strains as mechanism for 
adaptation. 

The second hypothesis investigated whether WEC 
could be used to detect adaptation in worms. We found 
that WEC over time was very noisy, even in the idealised 
situation of constant environmental conditions. There 
was evidence for adaptation by an increase in WEC over 
time, but in some parameter sets worms adapted by 
losing the favourable B allele, i.e. by reducing their abi- 
lity to parasitize sheep carrying the A allele but impro- 
ving their ability to survive on pasture. Thus, an increase 
in WEC over time did not necessarily imply that worms 
had improved their ability to parasitize sheep. The ob- 
served increase in WEC was also smaller than naively 
expected, i.e. WEC did not increase to the same level as 
in susceptible (x = 0) sheep, even when worms com- 
pletely fixed the favourable B allele. This occurred 
because the antagonistic pleiotropy of the worm B allele. 
Thus, although worms improved their ability to parasit- 
ise sheep by increasing the frequency of the B allele, they 
also reduced their survival on pasture, therefore lower- 
ing larval intake, worm burden and WEC of sheep. 

Our conclusions of course depend on the assumptions 
used in our model, i.e. the properties of the interacting 
loci in sheep and worms (Table 1), which represent our 
understanding of the processes involved in adaptation. 
For the source of genetic variation in worms, we make a 
distinction between segregating and de novo mutations, 
where the former has implicit restrictions on the effects 
of polymorphisms present in worm populations because 
of natural selection. We expect de novo mutations to 
add genetic variation each generation and, given the 
large effective population size of worms, this may be a 
substantial contribution. For example, Barton [28] ar- 
gued that it is possible to observe a specific nucleotide 
substitution in each generation of the worm population, 
assuming a population of 10 6 individuals (i.e. a mutation 
rate of 10" 8 /3, where 10" 8 is the mutation rate per 
nucleotide and there are possible three nucleotide 
substitutions). However the forces that act on de novo 
mutations are the same as we have modelled for segre- 
gating mutations, except that the allele frequency for de 
novo mutations is very low. Hence the selection differen- 
tial for the worm population in each sheep genotype and 
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the proportion of the worm population exposed to this 
selection pressure will determine if a new mutation per- 
sists, and the rate at which it will increase in frequency 
should it persist. However to facilitate differential adap- 
tation (i.e. adaptation specifically to sheep with A geno- 
type), these new mutations still require them to have 
favourable effects on worm fitness and increase in 
frequency in low WEC sheep but be unfavourable (or 
neutral) in sheep that are not selected. 

To complete our findings, several issues need to be 
further studied. The first limitation of our model is that 
we consider only a single interacting locus. Fishers 
infinitesimal model (a very large number of loci each 
with a negligible effect) and recent experimental evi- 
dence suggests that WEC in sheep is a highly polygenic 
trait [29]. It might be expected that presence of multiple 
loci affecting WEC would reduce the selection pressure 
acting on each worm locus, although this is yet to be 
rigorously explored. Secondly, a suitable method to test 
for adaptation, in theory and practice, might be to mea- 
sure WEC in a cross-factorial design with artificial in- 
fection e.g. [7]. This approach should limit the impact of 
antagonistic effects on the detection of adaptation 
because each sheep would be exposed to the same larval 
intake, but it does negate effects that accumulate during 
the season. Finally, the ability of worms to adapt to an- 
thelmintics should be comparable to adaptation under 
the framework we have established. It is well known that 
worms can adapt to anthelmintics e.g. [30] and at first 
sight this outcome differs from results presented here. 
However the selection pressures placed on worms by 
anthelmintic administration are much greater than those 
exerted by host genetic effects. These factors can be 
accounted for through the available parameters and are 
the subject of ongoing research. 

Conclusions 

Our model suggests that detectable adaptation of worms 
to sheep that are selected for low WEC is unlikely in the 
short term. This is because mutations in the worm that 
have properties suitable for adaptation, i.e. mutations 
that are favourable for survival in low WEC sheep but 
unfavourable in unselected sheep, may take many gener- 
ations to increase in frequency or are likely to be neutral 
(or near neutral) with respect to overall worm fitness in 
the current population. Even when we assumed ideal 
conditions and adaptation was observed in our model, 
the expected increase in WEC as a result of adaptation 
was small relative to the temporal fluctuations in the 
trait and also small relative to the reduction in WEC 
expressed by the sheep. The selection pressure applied 
to the worm population was a key factor for worm 
adaptation, with the selection differential between 



worm genotypes and the proportion of the worm 
population exposed to the selection pressure (i.e. the 
frequency of the low WEC allele in sheep) being key 
parameters. 

Additional file 



Additional file 1: Details of the host-parasite interaction model. 

This file contains all equations and parameters to reconstruct the full 
host-parasite interaction model. 
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